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ABSTRACT 


A  14.12-percent-thick,  natural-laminar-flow  airfoil,  the  S415,  intended  for  the  rotor  of 
a  helicopter  in  hover  has  been  designed  and  analyzed  theoretically  and  verified  experimen¬ 
tally  in  The  Pennsylvania  State  University  Low-Speed,  Low-Turbulence  Wind  Tunnel.  The 
two  primary  objectives  of  high  maximum  lift  and  low  profile  drag  have  been  achieved.  The 
constraint  on  the  pitching  moment  has  been  satisfied.  The  airfoil  exhibits  a  docile  stall.  Com¬ 
parisons  of  the  theoretical  and  experimental  results  generally  show  good  agreement. 


INTRODUCTION 


Almost  all  airfoils  in  use  on  rotorcraft  today  were  developed  under  the  assumption  that 
extensive  laminar  flow  is  not  likely  on  a  rotor.  (See  ref.  1,  for  example.)  For  the  present 
application,  however,  given  the  moderate  Reynolds  numbers  and  the  exploratory  nature  of  the 
current  effort,  the  achievement  of  laminar  flow  warrants  examination. 

The  airfoil  designed  under  the  present  effort  is  intended  for  the  rotor  of  a  helicopter  in 
hover.  The  airfoil  would  then  be  morphed  into  a  shape  more  suitable  for  cruise.  (See  ref.  2.) 
To  complement  the  design  effort,  an  investigation  was  conducted  in  The  Pennsylvania  State 
University  Low-Speed,  Low-Turbulence  Wind  Tunnel  (ref.  3)  to  obtain  the  basic,  low-speed, 
two-dimensional  aerodynamic  characteristics  of  the  airfoil.  The  results  have  been  compared 
with  predictions  from  the  method  of  references  4  and  5  and  from  the  method  of  reference  6. 


SYMBOLS 


Values  are  given  in  both  SI  and  U.S.  Customary  Units, 
tions  were  made  in  U.S.  Customary  Units. 


a 


pressure  coefficient, 


P/-Po 


Measurements  and  calcula- 


c 


c 


c 


airfoil  chord,  mm  (in.) 


section  chord-force  coefficient, 


cd 


section  profile-drag  coefficient, 
cnsina  +  eceosa 


Wake 


except  post  stall, 


C(j'  point  drag  coefficient  (ref.  7) 

c i  section  lift  coefficient,  cn/ cosa  -  cdtana 


c 


m 


section  pitching-moment  coefficient  about  quarter-chord  point, 
^Cp(5w°.25)d(i)4cp(5)d(5) 


section  nonnal-force  coefficient,  -  jjCpd(^- 


h  horizontal  width  in  wake  profile,  mm  (in.) 

M  free-stream  Mach  number 

p  static  pressure,  Pa  (lbf/ft  ) 

q  dynamic  pressure,  Pa  (lbf/ft  ) 

R  Reynolds  number  based  on  free-stream  conditions  and  airfoil  chord 

s  arc  length  along  airfoil  surface,  mm  (in.) 

t  airfoil  thickness,  mm  (in.) 

x  airfoil  abscissa,  mm  (in.) 

y  model  span  station,  y  =  0  at  midspan,  mm  (in.) 

z  airfoil  ordinate,  mm  (in.) 

a  angle  of  attack  relative  to  x-axis,  deg 

Subscripts: 

/  local  point  on  airfoil 

11  lower  limit  of  low-drag  range 

max  maximum 

min  minimum 

S  separation 

T  transition 

ul  upper  limit  of  low-drag  range 
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oo 


free-stream  conditions 


Abbreviations: 

L.  lower  surface 

NACA  National  Advisory  Committee  for  Aeronautics 

S.  boundary-layer  separation  location,  xs/c 

T.  boundary-layer  transition  location,  xT/c 

U.  upper  surface 


AIRFOIL  DESIGN 
OBJECTIVES  AND  CONSTRAINTS 

The  airfoil  design  specifications  are  contained  in  table  I.  Two  primary  objectives  are 
evident.  The  first  objective  is  to  achieve  a  maximum  lift  coefficient  of  1 .50  at  a  Mach  number 
of  0.50  and  a  Reynolds  number  of  5.00  x  106.  A  requirement  related  to  this  objective  is  that 
the  maximum  lift  coefficient  not  decrease  significantly  with  transition  fixed  near  the  leading 
edge  on  both  surfaces.  In  addition,  the  airfoil  should  exhibit  docile  stall  characteristics  at  a 
Mach  number  of  0.2  and  a  Reynolds  number  of  2.0  x  106,  which  is  within  the  operational 
range  of  the  Penn  State  Low-Speed,  Low-Turbulence  Wind  Tunnel.  The  second  objective  is 
to  obtain  low  profile-drag  coefficients  from  a  lift  coefficient  of  0.60  to  1.40  at  a  Mach  number 
of  0.50  and  a  Reynolds  number  of  5.00  x  106. 

One  major  constraint  was  placed  on  the  design  of  the  airfoil.  At  a  Mach  number  of 
0.40  and  a  Reynolds  number  of  4.00  x  106,  the  pitching-moment  coefficient  at  a  lift  coeffi¬ 
cient  of  1 .5  must  be  no  more  negative  than  -0.10. 
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PHILOSOPHY 


Given  the  above  objectives  and  constraint,  certain  characteristics  of  the  design  are 
apparent.  The  following  sketch  illustrates  a  drag  polar  that  meets  the  goals  for  this  design. 


cd 


Sketch  1 


The  desired  airfoil  shape  can  be  traced  to  the  pressure  distributions  that  occur  at  the  various 
points  in  sketch  1 .  Point  A  is  the  lower  limit  of  the  low-drag  range  of  lift  coefficients;  point  B, 
the  upper  limit.  The  profile-drag  coefficient  at  point  B  is  not  as  low  as  at  point  A,  unlike  the 
polars  of  many  laminar-flow  airfoils  where  the  drag  coefficient  within  the  laminar  bucket  is 
nearly  constant.  (See,  for  example,  ref.  8.)  This  characteristic  is  related  to  the  elimination  of 
significant  (i.e.,  drag-producing)  laminar  separation  bubbles  on  the  upper  surface  for  the 
design  range  of  Reynolds  numbers.  (See  ref.  9.)  The  drag  coefficient  increases  rapidly  out¬ 
side  the  low-drag,  lift-coefficient  range  because  boundary-layer  transition  moves  quickly 
toward  the  leading  edge  with  increasing  (or  decreasing)  lift  coefficient.  This  feature  results  in 
a  leading-edge  shape  that  produces  a  suction  peak  at  higher  lift  coefficients,  which  ensures 
that  transition  on  the  upper  surface  will  occur  very  near  the  leading  edge.  Thus,  the  maximum 
lift  coefficient,  point  C,  occurs  with  turbulent  flow  along  the  entire  upper  surface  and,  there¬ 
fore,  should  be  relatively  insensitive  to  roughness  at  the  leading  edge. 
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From  the  preceding  discussion,  the  pressure  distributions  along  the  polar  can  be 
deduced.  The  pressure  distribution  at  point  A  should  look  something  like  sketch  2. 


Sketch  2 


To  achieve  low  drag,  a  favorable  pressure  gradient  is  desirable  along  the  upper  surface  to 
about  30-percent  chord.  Aft  of  this  point,  a  short  region  having  a  shallow,  adverse  pressure 
gradient  (i.e.,  a  “transition  ramp”)  promotes  the  efficient  transition  from  laminar  to  turbulent 
flow  (ref.  10).  The  transition  ramp  is  followed  by  a  concave  pressure  recovery,  which  exhibits 
lower  drag  and  has  less  tendency  to  separate  than  the  corresponding  linear  or  convex  pressure 
recovery  (ref.  10).  The  specific  pressure  recovery  employed  represents  a  compromise 
between  maximum  lift,  drag,  pitching  moment,  stall  characteristics,  and  drag  divergence.  The 
steep,  adverse  pressure  gradient  aft  of  about  90-percent  chord  is  a  “separation  ramp,”  origi¬ 
nally  proposed  by  F.  X.  Wortmann,1  which  confines  turbulent  separation  to  a  small  region 
near  the  trailing  edge.  By  constraining  the  movement  of  the  separation  point  at  high  angles  of 
attack,  higher  lift  coefficients  can  be  achieved  with  little  drag  penalty.  (See  ref.  11.)  This  fea¬ 
ture  has  the  added  benefit  of  promoting  docile  stall  characteristics.  (See  ref.  12.) 

Along  the  lower  surface,  the  pressure  gradient  is  initially  adverse  and  then  favorable  to 
about  45-percent  chord.  Thus,  transition  is  imminent  over  the  forward  portion  of  the  lower 
surface.  (See  ref.  13.)  This  concept  allows  a  wide  low-drag  range  to  be  achieved  and 
increases  the  loading  in  the  leading-edge  region.  The  forward  loading  serves  to  balance,  with 
respect  to  the  pitching-moment  constraint,  the  aft  loading,  both  of  which  contribute  to  the 


'Director,  Institute  for  Aerodynamics  and  Gas  Dynamics,  University  of  Stuttgart,  Germany,  1974-1985. 
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achievement  of  a  high  maximum  lift  coefficient  and  low  profile-drag  coefficients.  This  region 
is  followed  by  a  transition  ramp  and  then  a  concave  pressure  recovery. 

The  amounts  of  pressure  recovery  on  the  upper  and  lower  surfaces  are  detennined  by 
the  width  of  the  low-drag,  lift-coefficient  range  and  the  pitching-moment  constraint. 

At  point  B,  the  pressure  distribution  should  look  like  sketch  3. 


Sketch  3 


No  suction  peak  exists  at  the  leading  edge.  Instead,  a  moderately  adverse  pressure  gradient 
extends  from  the  leading  edge  to  the  beginning  of  the  pressure  recovery. 


EXECUTION 

Given  the  pressure  distributions  previously  discussed,  the  design  of  the  airfoil  is 
reduced  to  the  inverse  problem  of  transforming  the  pressure  distributions  into  an  airfoil  shape. 
The  Eppler  Airfoil  Design  and  Analysis  Code  (refs.  4  and  5)  was  used  because  of  its  unique 
capability  for  multipoint  design  and  because  of  confidence  gained  during  the  design,  analysis, 
and  experimental  verification  of  many  other  airfoils.  (See  ref.  14,  for  example.) 
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The  airfoil  is  designated  the  S415.  The  airfoil  shape  and  coordinates  are  available 
from  Airfoils,  Incorporated.  The  airfoil  thickness  is  14.12-percent  chord. 


THEORETICAL  PROCEDURE 


The  theoretical  results  are  predicted  using  the  method  of  references  4  and  5 
(PROFIL07),  commonly  known  as  the  Eppler  code,  and  the  method  of  reference  6 
(MSES  3.0).  Critical  amplification  factors  of  11  and  9  were  specified  for  the  boundary-layer 
transition  computations  using  the  method  of  references  4  and  5  and  the  method  of  reference  6, 
respectively.  Because  the  maximum  lift  coefficient  computed  by  the  method  of  references  4 
and  5  is  not  always  realistic,  an  empirical  criterion  has  been  applied  to  the  computed  results. 
The  criterion  assumes  the  maximum  lift  coefficient  has  been  reached  if  the  drag  coefficient  of 
the  upper  surface  reaches  a  certain  value  that  is  a  function  of  the  Reynolds  number  and  the 
wind-tunnel  facility.  It  should  also  be  noted  that  the  compressibility  correction  (ref.  15)  incor¬ 
porated  in  the  method  of  references  4  and  5  is  invalid  if  the  local  flow  is  supersonic. 

Because  the  free-stream  Mach  number  for  all  wind-tunnel  test  conditions  did  not 
exceed  0.2,  the  flow  can  be  considered  essentially  incompressible  for  the  purpose  of  compar¬ 
ing  the  theoretical  and  experimental  results.  This  allows  the  (incompressible)  conformal- 
mapping  (design)  method  of  references  4  and  5  and  the  fast,  subcritical  flow  solver  of  the 
method  of  reference  6  to  be  used. 


EXPERIMENTAL  PROCEDURE 
WIND  TUNNEL 

The  Pennsylvania  State  University  Low-Speed,  Low-Turbulence  Wind  Tunnel  (ref.  3) 
is  a  closed-throat,  single-return,  atmospheric  tunnel  (fig.  1).  The  test  section  is  101.3  cm 
(39.9  in.)  high  by  147.6  cm  (58.1  in.)  wide  (fig.  2).  Electrically  actuated  turntables  provide 
positioning  and  attachment  for  the  two-dimensional  model.  The  turntables  are  flush  with  the 
top  and  bottom  tunnel  walls  and  rotate  with  the  model.  The  axis  of  rotation  coincided  approx¬ 
imately  with  midchord.  The  model  was  mounted  vertically  between  the  turntables  and  the 
gaps  between  the  model  and  the  turntables  were  sealed.  The  turbulence  intensity  in  the  test 
section  is  approximately  0.05  percent  at  46  m/s  (150  ft/s). 


MODEL 

The  aluminum,  wind-tunnel  model  was  fabricated  by  Advanced  Technologies,  Incor¬ 
porated,  Newport  News,  Virginia,  using  a  numerically  controlled  milling  machine.  The  model 
had  a  chord  of  533.38  mm  (20.999  in.)  and  a  span  of  107.95  cm  (42.50  in.)  and,  thus,  extended 
through  both  turntables.  Upper-  and  lower-surface  orifices  were  located  to  one  side  of  mid¬ 
span  at  the  staggered  positions  listed  in  table  II.  All  the  orifices  were  0.51  mm  (0.020  in.)  in 
diameter  with  their  axes  perpendicular  to  the  surface.  The  surfaces  of  the  model  were  sanded 
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to  ensure  an  aerodynamically  smooth  finish.  The  measured  model  contour  was  within 
0.13  mm  (0.005  in.)  of  the  prescribed  shape. 


WAKE-SURVEY  PROBE 

A  total-  and  static-pressure,  wake-survey  probe  (fig.  3)  was  mounted  from  the  top  tun¬ 
nel  wall  (fig.  2).  The  probe  was  positioned  50.8  cm  (20.0  in.)  from  the  ceiling  (i.e.,  about 
midspan)  and  automatically  aligned  with  the  wake-centerline  streamline.  A  traverse  mecha¬ 
nism  incrementally  positioned  the  probe  to  survey  the  wake.  The  increment  was  1 .27  mm 
(0.050  in.)  for  traverses  less  than  254.0  mm  (10.00  in.)  and  2.54  mm  (0.100  in.)  for  longer  tra¬ 
verses,  which  were  occasionally  required  near  the  maximum  lift  coefficient.  The  tip  of  the 
probe  was  located  0.6  chord  downstream  of  the  trailing  edge  of  the  model. 


INSTRUMENTATION 

Basic  tunnel  pressures  and  the  wake  pressures  were  measured  with  precision  transduc¬ 
ers.  Measurements  of  the  pressures  on  the  model  were  made  by  an  automatic  pressure¬ 
scanning  system  utilizing  precision  transducers.  Data  were  obtained  and  recorded  by  an  elec¬ 
tronic  data-acquisition  system. 


METHODS 

The  pressures  measured  on  the  model  were  reduced  to  standard  pressure  coefficients 
and  numerically  integrated  to  obtain  section  nonnal-force  and  chord-force  coefficients  and 
section  pitching-moment  coefficients  about  the  quarter-chord  point.  Section  profile-drag 
coefficients  were  computed  from  the  wake  total  and  static  pressures  by  the  method  of  refer¬ 
ence  7.  Wake  surveys  were  not  perfonned,  however,  at  most  post-stall  angles  of  attack,  in 
which  case,  the  profile-drag  coefficients  were  computed  from  the  nonnal-  and  chord-force 
coefficients. 

Standard,  low-speed,  wind-tunnel  boundary  corrections  (ref.  16)  have  been  applied  to 
the  data.  It  should  be  noted  that  the  model-chord-to-tunnel-height  ratio  is  0.36,  which  is 
approaching  the  limit  of  0.4  given  in  reference  16.  The  wake-survey-probe  total-pressure- 
tube  displacement  correction  (ref.  7)  has  been  taken  into  account. 


TESTS 

The  model  was  tested  at  Reynolds  numbers  based  on  airfoil  chord  of  1.0  x  1 06, 
1.5  x  106,  and  2.0  x  106  with  transition  free  (smooth)  and  with  transition  fixed  by  roughness 
near  the  leading  edge,  2-percent  chord  on  the  upper  surface  and  5-percent  chord  on  the  lower 
surface,  to  simulate  full-chord,  turbulent  flow.  The  grit  roughness  was  sized  using  the  method 
of  reference  17  and  sparsely  distributed  along  3-mm  (0.1 -in.)  wide  strips  applied  to  the  model 


with  lacquer.  (See  table  111(a).)  The  model  was  also  tested  with  a  roughness  equivalent  to 
NACA  standard  roughness  (ref.  8),  which  consisted  of  grit  roughness  having  a  nominal  size  of 
0.249  mm  (0.0098  in.)  applied  to  the  model  with  lacquer  and  sparsely  distributed  from  the 
leading  edge  to  an  arc  length  of  8-percent  chord  on  the  upper  and  lower  surfaces.  (See 
table  111(b).)  (The  grit  size  was  scaled  from  the  NACA  standard-roughness  grit  size  by  the 
ratio  of  the  model  chords  used  in  the  two  wind  tunnels:  533.4  mm  (21.00  in.)  in  the  present 
investigation  and  609.6  mm  (24.00  in.)  in  the  NACA  tests.) 

The  Mach  number  did  not  exceed  0.2  for  any  test  condition.  Thus,  the  test  Mach  num¬ 
bers  and  Reynolds  numbers  are  much  lower  than  the  operational  values  of  the  intended  appli¬ 
cation. 


Starting  from  1°,  which  is  within  the  low-drag  range  of  lift  coefficients,  the  angle  of 
attack  was  decreased  to  near  that  for  zero  lift.  The  angle  of  attack  was  then  increased  from  1° 
to  post-stall  values. 


DISCUSSION  OF  RESULTS 
THEORETICAL  RESULTS 
Pressure  Distributions 

The  pressure  distributions  predicted  using  the  method  of  reference  6  (MSES  3.0)  at 
various  angles  of  attack  at  a  Mach  number  of  0.50  and  a  Reynolds  number  of  5.00  x  106  are 
shown  in  figure  4. 


Section  Characteristics 

The  section  characteristics  predicted  using  the  method  of  references  4  and  5 
(PROFIL07)  and  the  method  of  reference  6  (MSES  3.0)  at  a  Mach  number  of  0.50  and  a 
Reynolds  number  of  5.00  x  106  with  transition  free  and  with  transition  fixed  are  shown  in  fig¬ 
ure  5.  (It  should  be  remembered  that  the  compressibility  correction  incorporated  in  the 
method  of  refs.  4  and  5  is  invalid  if  the  local  flow  is  supersonic  and,  accordingly,  only  sub¬ 
sonic  results  are  shown.)  Based  on  the  predictions,  all  the  design  objectives  and  the  design 
constraint  have  essentially  been  met. 


EXPERIMENTAL  RESULTS 
Pressure  Distributions 

The  pressure  distributions  at  various  angles  of  attack  for  a  Reynolds  number  of 
1.5  x  106  and  a  Mach  number  of  0. 1  with  transition  free  are  shown  in  figure  6.  At  an  angle  of 
attack  of -5. 10°  (fig.  6(a)),  a  short  laminar  separation  bubble  is  discernible  on  the  upper  sur- 
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face  around  55-percent  chord  and  a  long  laminar  separation  bubble  is  evident  on  the  lower 
surface  aft  of  the  leading  edge.  As  the  angle  of  attack  is  increased,  the  short  bubble  on  the 
upper  surface  moves  forward  and  the  long  bubble  on  the  lower  surface  decreases  in  length, 
until  it  has  disappeared  at  an  angle  of  attack  of -2.02°  (fig.  6(a)).  At  an  angle  of  attack  of 
0.03°  (fig.  6(b)),  a  short  laminar  separation  bubble  is  evident  on  the  lower  surface  around 
65-percent  chord,  the  laminar  flow  having  survived  the  pressure  peak  near  the  leading  edge. 
As  the  angle  of  attack  is  increased  further,  the  bubble  on  the  upper  surface  continues  to  move 
forward,  whereas  the  bubble  on  the  lower  surface  migrates  slowly  aft  (figs.  6(b)-6(d)).  At  an 
angle  of  attack  of  8.23°  (fig.  6(d)),  turbulent,  trailing-edge  separation  is  evident  on  the  upper 
surface.  The  amount  of  separation  increases  with  increasing  angle  of  attack  (fig.  6(d)).  The 
maximum  lift  coefficient  occurs  at  an  angle  of  attack  just  beyond  9.24°  (fig.  6(d)).  As  the 
angle  of  attack  is  increased  even  further,  the  separation  point  continues  to  move  forward  and 
the  leading-edge  peak  continues  to  climb  (fig.  6(e)). 


Section  Characteristics 

The  section  characteristics  with  transition  free,  with  transition  fixed,  and  with  scaled, 
NACA  standard  roughness,  denoted  “rough,”  are  shown  in  figure  7  and  tabulated  in  the 
appendix.  For  a  Reynolds  number  of  2.0  x  106  and  a  Mach  number  of  0.2  with  transition  free 
(fig.  7(c)),  the  maximum  lift  coefficient  is  1.37  and  the  stall  characteristics  are  relatively  doc¬ 
ile.  The  lower  limit  of  the  low-drag,  lift-coefficient  range  is  0.50;  the  upper  limit,  1.28.  For  a 
Reynolds  number  of  1.5  x  106  and  a  Mach  number  of  0.1  with  transition  free  (fig.  7(b)),  the 
zero-lift  pitching-moment  coefficient  is  -0.09. 

The  effects  of  Reynolds  number  on  the  section  characteristics  are  summarized  in  fig¬ 
ure  8.  In  general,  the  lift-curve  slope,  the  maximum  lift  coefficient,  the  lower  limit  of  the  low- 
drag  range,  and  the  magnitude  of  the  pitching-moment  coefficients  increase  with  increasing 
Reynolds  number  and  the  stall  characteristics  become  less  docile.  The  profile-drag  coeffi¬ 
cients  and  the  upper  limit  of  the  low-drag  range  decrease  with  increasing  Reynolds  number. 

The  effect  of  fixing  transition  on  the  section  characteristics  is  shown  in  figure  7.  In 
general,  the  lift-curve  slope,  the  maximum  lift  coefficient,  and  the  magnitude  of  the  pitching- 
moment  coefficients  decrease  with  transition  fixed.  These  results  are  primarily  a  consequence 
of  the  boundary-layer  displacement  effect,  which  decambers  the  airfoil  because  the  displace¬ 
ment  thickness  is  greater  with  transition  fixed  than  with  transition  free.  In  addition,  the  maxi¬ 
mum  lift  coefficient  decreases  with  transition  fixed  because  the  roughness  induces  earlier 
trailing-edge  separation.  The  reduction  in  maximum  lift  coefficient  averages  4  percent  over 
the  test  Reynolds  number  range.  The  drag  coefficients  are,  of  course,  adversely  affected  by 
the  roughness. 

It  should  be  noted  that,  for  many  test  conditions,  the  Reynolds  number  based  on  local 
velocity  and  boundary-layer  displacement  thickness  at  the  roughness  locations  is  too  low  to 
support  turbulent  flow.  (See  ref.  18.)  Accordingly,  to  force  transition,  the  roughness  must  be 
so  large  that  it  increases  the  displacement  thickness,  which  abnormally  decreases  the  lift  coef¬ 
ficient  and  the  magnitude  of  the  pitching-moment  coefficient  and  increases  the  drag  coefifi- 


10 


cient.  Conversely,  at  low  lift  coefficients,  the  roughness  on  the  upper  surface,  which  is  sized 
for  high  lift  coefficients,  is  too  small  to  force  transition,  resulting  in  incorrectly  low  drag  coef¬ 
ficients. 

The  effect  of  the  scaled,  NACA  standard  roughness  on  the  section  characteristics  is 
shown  in  figure  7.  The  effect  is  more  severe  than  that  of  fixing  transition.  The  reduction  in 
maximum  lift  coefficient  is  larger,  averaging  13  percent  over  the  test  Reynolds  number  range. 
It  should  be  remembered  that  the  effect  of  roughness  is  proportional  to  the  ratio  of  the  rough¬ 
ness  height  to  the  boundary-layer  thickness.  Because  the  height  of  the  scaled,  NACA  stan¬ 
dard  roughness  and  the  airfoil  chord  are  constant,  the  effect  of  this  roughness  typically 
increases  with  increasing  Reynolds  number  (because  increasing  Reynolds  number  results  in 
decreasing  boundary-layer  thickness),  although  the  effect  on  the  maximum  lift  coefficient  is 
almost  constant  over  the  test  Reynolds  number  range. 

The  variations  of  maximum  lift  coefficient  and  minimum  profile-drag  coefficient  with 
Reynolds  number  are  shown  in  figures  9  and  10,  respectively.  In  general,  the  maximum  lift 
coefficient  increases  with  increasing  Reynolds  number,  whereas  the  minimum  profile-drag 
coefficient  decreases,  which  are  typical  trends  for  most  airfoils. 


COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS 

Pressure  Distributions 

The  comparison  of  the  theoretical  and  experimental  pressure  distributions  at  various 
angles  of  attack  is  shown  in  figure  11.  It  should  be  noted  that  the  pressure  distributions  pre¬ 
dicted  using  the  method  of  references  4  and  5  (PROFIL07)  are  inviscid  and  incompressible, 
whereas  the  pressure  distributions  predicted  using  the  method  of  reference  6  (MSES  3.0)  as 
well  as  the  experimental  pressure  distributions  were  obtained  for  a  Reynolds  number  of 
1.5  x  106  and  a  Mach  number  of  0. 1  with  transition  free.  It  should  also  be  noted  that  the  theo¬ 
retical  lift  coefficient  from  the  method  of  references  4  and  5  is  calculated  from  the  lift-curve 
slope  and  the  angle  of  attack  relative  to  the  zero-lift  line,  whereas  the  lift  coefficient  from  the 
method  of  reference  6  and  from  the  experiment  is  derived  from  the  integrated  pressure  distri¬ 
bution.  (See  refs.  4-7.)  Thus,  at  a  given  lift  coefficient,  the  pressure  distribution  predicted 
using  the  method  of  references  4  and  5  does  not  necessarily  have  the  same  area  as  the  mea¬ 
sured  pressure  distribution. 

With  respect  to  the  method  of  references  4  and  5,  at  lift  coefficients  of  0.48  and  0.93 
(figs.  11(a)  and  11(b),  respectively),  the  pressure  coefficients  and  the  pressure  gradients  agree 
reasonably  well,  except  in  the  vicinity  of  the  laminar  separation  bubbles  and  near  the  trailing 
edge.  The  latter  disparity  is  probably  the  result  of  the  boundary-layer  displacement  effect.  At 
a  lift  coefficient  of  1.33  (fig.  11(c)),  which  is  just  below  the  measured  maximum  lift  coeffi¬ 
cient,  the  agreement  is  poor  because  the  effect  of  the  upper-surface,  trailing-edge  separation 
on  the  pressure  distribution  is  not  modelled  in  the  method  of  references  4  and  5. 
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With  respect  to  the  method  of  reference  6,  at  a  lift  coefficient  of  0.48  (fig.  11(a)),  the 
pressure  coefficients  and  the  pressure  gradients  agree  reasonably  well,  except  near  the  trailing 
edge  on  the  upper  surface.  The  location  of  the  upper-surface  laminar  separation  bubble  is  pre¬ 
dicted  well,  but  that  of  the  lower-surface  bubble  is  forward  of  the  measured  location.  At  a  lift 
coefficient  of  0.93  (fig.  1 1(b)),  the  pressure  distributions  again  agree  reasonably  well,  and  the 
locations  of  the  bubbles  are  predicted  well.  At  a  lift  coefficient  of  1.33  (fig.  11(c)),  the  agree¬ 
ment  is  poor  because  the  significant,  upper-surface,  trailing-edge  separation  is  not  predicted 
by  the  method  of  reference  6. 


Section  Characteristics 

The  comparison  of  the  theoretical  and  experimental  section  characteristics  with  transi¬ 
tion  free  is  shown  in  figure  12.  The  maximum  lift  coefficient  is  estimated  within  approxi¬ 
mately  1  percent  using  the  empirical  criterion  previously  discussed  applied  to  the  results  from 
the  method  of  references  4  and  5  (PROFIL07).  The  method  of  reference  6  (MSES  3.0)  over¬ 
predicts  the  maximum  lift  coefficient  by  an  average  of  13  percent.  The  method  of  reference  6 
underpredicts  the  profile-drag  coefficients.  The  method  of  references  4  and  5  underpredicts 
and  the  method  of  reference  6  overpredicts  the  upper  limit  of  the  low-drag,  lift-coefficient 
range.  The  method  of  references  4  and  5  overpredicts  the  magnitude  of  the  pitching-moment 
coefficients. 

The  comparison  of  the  theoretical  and  experimental  section  characteristics  with  transi¬ 
tion  fixed  is  shown  in  figure  13.  In  general,  the  predicted  characteristics  show  similar  tenden¬ 
cies  as  with  transition  free,  although  the  general  agreement  is  poorer,  probably  because  of  the 
abnonnalities  introduced  by  the  roughness,  as  discussed  previously.  In  particular,  the  method 
of  reference  6  overpredicts  the  maximum  lift  coefficient  by  an  average  of  18  percent. 


CONCLUDING  REMARKS 


A  14.12-percent-thick,  natural-laminar- flow  airfoil,  the  S415,  intended  for  the  rotor  of 
a  helicopter  in  hover  has  been  designed  and  analyzed  theoretically  and  verified  experimen¬ 
tally  in  The  Pennsylvania  State  University  Low-Speed,  Low-Turbulence  Wind  Tunnel.  The 
two  primary  objectives  of  a  high  maximum  lift  coefficient  and  low  profile-drag  coefficients 
have  been  achieved.  The  constraint  on  the  pitching-moment  coefficient  has  been  satisfied. 
The  airfoil  exhibits  docile  stall  characteristics.  Comparisons  of  the  theoretical  and  experimen¬ 
tal  results  generally  show  good  agreement. 
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TABLE  I.-  AIRFOIL  DESIGN  SPECIFICATIONS 


Parameter 

Objective/ 

Constraint 

Mach 

Number 

M 

Reynolds 

Number 

R 

Priority 

Minimum  lift  coefficient 

G,min 

0.40 

Medium 

Maximum  lift  coefficient 

G,max 

1.50 

0.50 

5.00  x  106 

High 

Lower  limit  of  low-drag, 
lift-coefficient  range  C/jj 

0.60 

Medium 

Upper  limit  of  low-drag, 
lift-coefficient  range  C/Uj 

1.40 

High 

Pitching-moment 
coefficient  cm  at  C/=  1.5 

>-0.10 

0.40 

4.00  x  106 

Low 

Thickness  t/c 

— 

— 

Other: 

Maximum  lift  coefficient  C/max  relatively  independent  of  leading-edge  roughness 
Docile  stall  characteristics  at  M  =  0.2  and  R  =  2.0xl06  (i.e.,  verifiable  in  tunnel) 


15 


TABLE  II.-  MODEL  ORIFICE  LOCATIONS 


[c  =  533.38  mm  (20.999  in.)] 


Upper  Surface 


Lower  Surface 


x/c 

y,  mm  (in.) 

x/c 

y,  mm  (in.) 

0.00160 

-129.31  (-5.091) 

0.00000 

-149.16  (-5.873) 

.00814 

-128.44  (-5.057) 

.00319 

-148.14  (-5.832) 

.01876 

-127.47  (-5.019) 

.01408 

-146.87  (-5.782) 

.03339 

-125.86  (-4.955) 

.03185 

-145.14  (-5.714) 

.05237 

-124.29  (-4.893) 

.05515 

-142.86  (-5.624) 

.07463 

-121.78  (-4.795) 

.08426 

-140.01  (-5.512) 

.10121 

-119.35  (-4.699) 

.11951 

-136.88  (-5.389) 

.13064 

-116.40  (-4.582) 

.15934 

-132.87  (-5.231) 

.16370 

-112.93  (-4.446) 

.20402 

-128.43  (-5.056) 

.19875 

-109.47  (-4.310) 

.25191 

-123.96  (-4.880) 

.23670 

-106.06  (-4.176) 

.30051 

-119.28  (-4.696) 

.27655 

-102.26  (-4.026) 

.34825 

-115.02  (-4.528) 

.31781 

-98.07  (-3.861) 

.39710 

-109.98  (-4.330) 

.36033 

-94.12  (-3.705) 

.44586 

-105.13  (-4.139) 

.40492 

-89.77  (-3.534) 

.49582 

-100.36  (-3.951) 

.45107 

-85.05  (-3.349) 

.54463 

-95.54  (-3.761) 

.49838 

-80.48  (-3.169) 

.59492 

-90.52  (-3.564) 

.54661 

-76.02  (-2.993) 

.64423 

-85.74  (-3.376) 

.59565 

-78.60  (-3.094) 

.69319 

-81.02  (-3.190) 

.64354 

-81.49  (-3.208) 

.74148 

-76.15  (-2.998) 

.69157 

-84.54  (-3.328) 

.78943 

-80.98  (-3.188) 

.73792 

-87.24  (-3.435) 

.83365 

-85.50  (-3.366) 

.78196 

-89.99  (-3.543) 

.87502 

-89.92  (-3.540) 

.82338 

-92.42  (-3.638) 

.91162 

-93.86  (-3.695) 

.86120 

-94.59  (-3.724) 

.94305 

-97.11  (-3.823) 

.89419 

-96.19  (-3.787) 

.96748 

-99.59  (-3.921) 

.92556 

-98.30  (-3.870) 

.98520 

-101.35  (-3.990) 

.95187 

-99.90  (-3.933) 

.99667 

-103.88  (-4.090) 

.97271 

-101.13  (-3.981) 

.98764 

-101.81  (-4.008) 

.99748 

-102.30  (-4.027) 

1.00000 

-96.63  (-3.804) 
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TABLE  III.-  ROUGHNESS  LOCATIONS  AND  SIZES 


(a)  Transition  Fixed 


R 

Upper  surface 

Lower  surface 

x/c 

Grit 

number 

Nominal  size, 
mm  (in.) 

x/c 

Grit 

number 

Nominal  size, 
mm  (in.) 

1.0  x  106 

0.02 

90 

0.178  (0.0070) 

0.05 

46 

0.419(0.0165) 

1.5  x  106 

100 

0.150  (0.0059) 

60 

0.297  (0.0117) 

2.0  x  106 

150 

0.104  (0.0041) 

70 

0.249  (0.0098) 

(b)  Scaled,  NACA  Standard  Roughness 


R 

Upper  surface 

Lower  surface 

s/c 

Grit 

number 

Nominal  size, 
mm  (in.) 

s/c 

Grit 

number 

Nominal  size, 
mm  (in.) 

1.0  x  106 

0  to 
0.08 

70 

0.249  (0.0098) 

0  to 
0.08 

70 

0.249  (0.0098) 

1.5  x  106 

2.0  x  106 
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motor 


Figure  1.-  The  Pennsylvania  State  University  Low-Speed,  Low-Turbulence  Wind  Tunnel. 


Figure  2.-  S415  airfoil  model  and  wake-survey  probe  mounted  in  test  section. 


Static-pressure  connection 
Total-pressure  connection 


Figure  3.-  Wake-survey  probe. 


Figure  4.-  Concluded. 


PR0FIL07 
MSES  3.0 


Figure  5.-  Concluded. 
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-4.08° 

-3.05° 

-2.02° 


24 

0 


0.5  x/c 


(a)  a  =  -5.10°,  -4.08°,  -3.05°,  and  -2.02°. 

Figure  6.- Experimental  pressure  distributions  for  R=1.5x  106  and  M  =  0.1 

free. 
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1 

1 


with  transition 


(d)  a  =  7.21°,  8.23°,  9.24°,  and  10.24°. 


Figure  6.-  Continued. 


Figure  6.-  Concluded. 


- ©  Transition  free 

- Q  Transition  fixed 


(a)  R  =  1.0  x  106  and  M  =  0.1. 

Figure  7.-  Experimental  section  characteristics  with  transition  free,  with  transition  fixed,  and  rough. 


Transition  free 


(b)  R=  1.5  x  106  and  M  =  0.1. 


Figure  7.-  Continued. 


(b)  Transition  fixed. 


Figure  8.-  Continued. 
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Figure  9.-  Variation  of  experimental  maximum  lift  coefficient  with  Reynolds  number. 
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Figure  10.-  Variation  of  experimental  minimum  profile-drag  coefficient  with  Reynolds  num¬ 
ber. 
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-  Theory  (PR0FIL07) 

- Theory  (MSES  3.0) 

© - ©  Experiment 
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'  ©  —  ©  2^ 


1®-©-©  0 


a  relative  to  the  x-ax'is 


(a)  c/  =  0.48. 


Figure  11.-  Comparison  of  theoretical  and  experimental  pressure  distributions  for 
R  =  1.5  x  106  and  M  =  0.1  with  transition  free. 


-  Theory  (PR0FIL07) 

- Theory  (MSES  3.0) 

© - ©  Experiment 


a  relative  to  the  x-axis 


(b)  C/  =  0.93. 
Figure  11.-  Continued. 
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(c)  c/=  1.33. 
Figure  11.-  Concluded. 
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(b)  R=  1.5  x  10' 


Figure  12.- 


and  M  =  0.1. 
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(a)  R=  1.0  x  106  and  M  =  0.1. 


Figure  13.-  Comparison  of  theoretical  and  experimental  section  characteristics  with  transition  fixed. 


APPENDIX 

EXPERIMENTAL  SECTION  CHARACTERISTICS 
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R  =  1.00  x  106,  M  =  0.08,  transition  free 


a,  deg 

c  / 

-4.072 

0.0482 

-3.047 

.1506 

-2.020 

.2547 

-.994 

.3634 

-.484 

.4148 

-.228 

.4458 

.027 

.4688 

.282 

.4933 

.538 

.5192 

1.050 

.5757 

2.076 

.6876 

3.100 

.7950 

4.125 

.9082 

5.150 

1.0159 

6.176 

1.1282 

7.199 

1.2262 

7.455 

1.2566 

7.711 

1.2787 

7.967 

1.2985 

8.222 

1.3125 

9.237 

1.3218 

9.742 

1.3177 

1.245 

1.3056 

11.227 

1.2445 

12.220 

1.1840 

13.217 

1.1627 

Cd 


0.026784 

-0.08971 

.018508 

-.08828 

.010751 

-.08468 

.009222 

-.08450 

.008906 

-.08577 

.008216 

-.08743 

.007531 

-.08783 

.007475 

-.08812 

.007607 

-.08855 

.007712 

-.08930 

.007909 

-.09030 

.008368 

-.09114 

.008532 

-.09186 

.009123 

-.09210 

.009564 

-.09225 

.009719 

-.09266 

.009626 

-.09317 

.009972 

-.09263 

.009992 

-.09123 

.010509 

-.08920 

.015533 

-.07582 

.019722 

-.06904 

.025156 

-.06352 

.035996 

-.06816 

.043106 

-.06069 

.055247 

-.05940 
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R  =  1.00  x  106,  M  =  0.08,  transition  fixed 


a,  deg 

c  / 

-4.072 

0.0528 

-3.047 

.1540 

-2.019 

.2592 

-.994 

.3657 

.029 

.4765 

1.052 

.5749 

2.076 

.6728 

3.099 

.7672 

4.123 

.8602 

5.144 

.9436 

6.167 

1.0352 

7.188 

1.1110 

8.208 

1.1833 

8.716 

1.2102 

9.224 

1.2356 

9.731 

1.2520 

1.235 

1.2722 

1.732 

1.2178 

11.210 

1.1958 

12.207 

1.1292 

13.209 

1.1329 

Cd 


0.026468 

-0.08995 

.018252 

-.08843 

.010675 

-.08513 

.009408 

-.08518 

.008968 

-.08716 

.009035 

-.08727 

.009741 

-.08620 

.010232 

-.08479 

.012384 

-.08298 

.014678 

-.08059 

.015765 

-.07911 

.017302 

-.07595 

.018755 

-.07252 

.020696 

-.07027 

.022488 

-.06796 

.026197 

-.06572 

.023339 

-.06586 

.109652 

-.05614 

.126229 

-.07527 

.159061 

-.06086 

.191461 

-.06061 

50 


R  =  1.00  x  106,  M  =  0.08,  rough 


a,  deg 

c  / 

-4.070 

0.0488 

-3.048 

.1331 

-2.021 

.2309 

-.996 

.3322 

.026 

.4325 

1.049 

.5287 

2.072 

.6281 

3.095 

.7292 

4.118 

.8238 

5.141 

.9149 

6.162 

.9960 

7.181 

1.0646 

8.197 

1.1075 

8.701 

1.1173 

9.200 

1.1293 

9.695 

1.1327 

1.191 

1.1186 

11.186 

1.0635 

12.196 

1.0893 

13.197 

1.0911 

Cd 


0.024084 

-0.08768 

.017182 

-.08460 

.013585 

-.08046 

.011917 

-.07890 

.011350 

-.08000 

.011819 

-.07981 

.012360 

-.07952 

.012832 

-.07961 

.013584 

-.07920 

.014328 

-.07691 

.016132 

-.07464 

.017912 

-.07113 

.021345 

-.06542 

.022813 

-.06412 

.022682 

-.06805 

.024964 

-.07420 

.028474 

-.07525 

.050962 

-.06749 

.059921 

-.06389 

.072355 

-.06367 
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R  =  1.50  x  106,  M  =  0.12,  transition  free 


a,  deg 

c  / 

-5.098 

-0.0618 

-4.075 

.0430 

-3.048 

.1518 

-2.019 

.2599 

-.994 

.3728 

-.226 

.4536 

.028 

.4832 

.284 

.5107 

.540 

.5394 

1.052 

.5907 

2.078 

.7028 

3.103 

.8177 

4.129 

.9287 

5.154 

1.0411 

6.179 

1.1516 

6.440 

1.1954 

6.692 

1.2051 

6.949 

1.2319 

7.210 

1.2742 

7.465 

1.2923 

8.226 

1.3127 

9.240 

1.3280 

9.746 

1.3315 

1.245 

1.3098 

11.221 

1.2529 

12.216 

1.1772 

13.217 

1.1731 

Cd 

cm 

0.037471 

-0.08927 

.027653 

-.09100 

.018100 

-.08958 

.009470 

-.08540 

.008470 

-.08667 

.008111 

-.08788 

.006646 

-.09059 

.006787 

-.09093 

.006769 

-.09115 

.006940 

-.09124 

.007191 

-.09177 

.007501 

-.09320 

.007771 

-.09407 

.008214 

-.09494 

.008233 

-.09543 

.008366 

-.09524 

.008646 

-.09514 

.008784 

-.09459 

.008888 

-.09363 

.009420 

-.09229 

.012033 

-.08652 

.017959 

-.07519 

.022211 

-.07042 

.041694 

-.06615 

.041132 

-.07806 

.043370 

-.06449 

.054974 

-.06248 

52 


R=  1.50  x  106,  M  =  0.13,  transition  fixed 


a,  deg 

c  / 

-5.097 

-0.0580 

-4.074 

.0493 

-3.047 

.1565 

-2.018 

.2648 

-.993 

.3715 

.030 

.4776 

1.054 

.5719 

2.077 

.6659 

3.101 

.7630 

4.125 

.8620 

5.148 

.9561 

6.171 

1.0463 

7.194 

1.1308 

8.214 

1.2049 

9.231 

1.2621 

9.734 

1.2688 

1.221 

1.2524 

11.212 

1.2100 

12.213 

1.1731 

13.215 

1.1644 

Cd 

cm 

0.037115 

-0.08926 

.027177 

-.09114 

.017795 

-.08953 

.009377 

-.08562 

.008756 

-.08554 

.010057 

-.08684 

.010765 

-.08492 

.011320 

-.08391 

.012030 

-.08266 

.012854 

-.08159 

.014054 

-.08001 

.014929 

-.07800 

.016704 

-.07564 

.018722 

-.07257 

.023754 

-.06914 

.038114 

-.06771 

.029319 

-.07760 

.044823 

-.07730 

.047744 

-.06700 

.056824 

-.06313 
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R  =  1.50  x  106,  M  =  0.13,  rough 


a,  deg 

c  / 

-5.097 

-0.0690 

-4.073 

.0295 

-3.047 

.1299 

-2.021 

.2301 

-.997 

.3358 

.027 

.4422 

1.051 

.5401 

2.075 

.6429 

3.099 

.7445 

4.123 

.8409 

5.146 

.9328 

6.168 

1.0148 

7.187 

1.0839 

8.202 

1.1306 

9.202 

1.1570 

9.697 

1.1476 

1.188 

1.1287 

11.185 

1.0675 

12.195 

1.0819 

13.198 

1.1008 

Cd 

cm 

0.035415 

-0.08669 

.025725 

-.08484 

.018999 

-.08215 

.012075 

-.08009 

.011274 

-.08050 

.011104 

-.08109 

.011518 

-.08073 

.011999 

-.08065 

.012738 

-.08005 

.013468 

-.07865 

.015002 

-.07709 

.016250 

-.07428 

.018606 

-.07086 

.021597 

-.06648 

.022713 

-.07360 

.025234 

-.07727 

.047662 

-.08210 

.050181 

-.07038 

.056831 

-.06423 

.071990 

-.06618 

54 


R  =  2.00  x  106,  M  =  0.17,  transition  free 


a,  deg 

c  / 

-4.078 

0.0414 

-3.049 

.1575 

-2.020 

.2681 

-.993 

.3867 

-.223 

.4738 

.032 

.5048 

.288 

.5333 

1.054 

.6026 

2.080 

.7214 

3.106 

.8366 

4.133 

.9504 

5.159 

1.0691 

6.185 

1.1817 

6.698 

1.2351 

7.211 

1.2824 

7.467 

1.3019 

8.231 

1.3278 

9.247 

1.3594 

9.753 

1.3654 

1.250 

1.3353 

11.222 

1.2454 

12.217 

1.1816 

13.220 

1.1897 

cd 

cm 

0.028329 

-0.09266 

.018045 

-.09098 

.008795 

-.08789 

.008093 

-.08852 

.007771 

-.08969 

.006335 

-.09189 

.006421 

-.09266 

.006472 

-.09319 

.006890 

-.09530 

.007314 

-.09570 

.007468 

-.09624 

.007743 

-.09742 

.008027 

-.09805 

.008238 

-.09765 

.008574 

-.09623 

.009882 

-.09494 

.013974 

-.08660 

.018968 

-.07755 

.022104 

-.07305 

.031283 

-.06867 

.038509 

-.07700 

.044925 

-.06617 

.057890 

-.06571 

55 


R  =  1.99  x  106,  M  =  0.17,  transition  fixed 


a,  deg 

c  / 

-3.049 

0.1558 

-2.020 

.2660 

-.994 

.3781 

.030 

.4704 

1.053 

.5631 

2.077 

.6658 

3.102 

.7711 

4.128 

.8757 

5.152 

.9733 

6.176 

1.0677 

7.199 

1.1563 

8.220 

1.2331 

8.737 

1.2870 

9.240 

1.2954 

9.744 

1.3030 

1.233 

1.2935 

11.216 

1.2271 

12.216 

1.1741 

13.219 

1.1880 

Cd 

cm 

0.018201 

-0.09108 

.008846 

-.08789 

.008986 

-.08760 

.010368 

-.08504 

.010774 

-.08372 

.011394 

-.08422 

.012112 

-.08369 

.012710 

-.08299 

.014005 

-.08185 

.014881 

-.08009 

.016784 

-.07775 

.018639 

-.07441 

.022409 

-.07077 

.023975 

-.06971 

.027872 

-.06682 

.030545 

-.07621 

.045895 

-.07927 

.046306 

-.06534 

.062936 

-.06733 

56 


R  =  2.00  x  106,  M  =  0.18,  rough 


a,  deg 

c  / 

-4.075 

0.0258 

-3.048 

.1343 

-2.022 

.2377 

-.996 

.3514 

.031 

.4659 

1.054 

.5613 

2.077 

.6628 

3.103 

.7694 

4.128 

.8708 

5.151 

.9632 

6.175 

1.0499 

7.195 

1.1198 

8.211 

1.1681 

9.206 

1.1869 

9.695 

1.1591 

1.187 

1.1127 

11.188 

1.0682 

12.198 

1.0942 

13.201 

1.1102 

cd 

cm 

0.028500 

-0.08610 

.019130 

-.08422 

.011424 

-.08288 

.010999 

-.08322 

.010943 

-.08353 

.011416 

-.08325 

.011947 

-.08357 

.012714 

-.08306 

.013518 

-.08201 

.015081 

-.08013 

.016448 

-.07730 

.019168 

-.07348 

.022659 

-.06895 

.023375 

-.07838 

.038048 

-.08493 

.047673 

-.08097 

.051510 

-.06974 

.059991 

-.06574 

.073498 

-.06685 
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